Portable impedance-sensing device for microorganism characterization in the field

A variety of biosensors have been proposed to quickly detect and measure the properties of individual microorganisms among heterogeneous populations, but challenges related to cost, portability, stability, sensitivity, and power consumption limit their applicability. This study proposes a portable microfluidic device based on impedance flow-cytometry and electrical impedance spectroscopy that can detect and quantify the size of microparticles larger than 45 µm, such as algae and microplastics. The system is low cost ($300), portable (5 cm \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}× 5 cm), low-power (1.2 W), and easily fabricated utilizing a 3D-printer and industrial printed circuit board technology. The main novelty we demonstrate is the use of square wave excitation signal for impedance measurements with quadrature phase-sensitive detectors. A linked algorithm removes the errors associated to higher order harmonics. After validating the performance of the device for complex impedance models, we used it to detect and differentiate between polyethylene microbeads of sizes between 63 and 83 µm, and buccal cells between 45 and 70 µm. A precision of 3% is reported for the measured impedance and a minimum size of 45 µm is reported for the particle characterization.


Principles of impedance-flow cytometry
The characterization of microorganisms can be performed using their impedance spectrum, which is a function of the resistivity, dielectric constant, and geometry of the substance under test (SUT). The resistivity and dielectric constant depends on the mobility and quantity of charge carriers in the material 41 . The impedance is defined from the complex Ohm's law based on the ratio of a voltage signal to a current signal 37 , as shown in Eq. 1, where Z is the impedance, V is the applied (or measured) voltage, I is the measured (or applied) current, φ is the phase difference between the voltage and current, ω is the angular excitation frequency, and j is the imaginary unit value. EIS measurement consists in injecting an AC sinusoidal waveform of a known voltage or current to the SUT and measuring its respective output current or voltage response, for a certain number of frequency.
The permittivity of the medium, membrane, and cytoplasm are fundamental properties of cells, making them suitable for microorganism characterization. The simplest model for a biological cell with a lipid layer plasma membrane is the single-shell model 28,[42][43][44][45][46] modelled by a homogeneous phase cytoplasm and an insulating thin shell. Bacteria, yeast, and plant cells are examples of microorganisms that possess a cell wall outside of the plasma membrane. This addition modifies the single shell model to a double shell model. An in-depth cells analysis and modeling strategy is available in (Asami, 2002) 42 .
IFC 41,44,46,47 is the technique offering the best results in the literature for whole cell characterization so far. It is a label-free non-invasive impedimetric measurement technique based on the Coulter counter 48 to measure the volume displacement of particles flowing in a fluid. The particles are detected using IFC by monitoring the impedance changes observed each time a particle passes over electrodes in a narrow channel. This is the case since the particles and fluid have different impedances. Figure 1 describes the principles of IFC. Electrodes are positioned on the walls of a microchannel and the impedance of a liquid flowing within is measured for several excitation frequency. A pulsed waveform resembling the one shown at the bottom of Fig. 1 is thus retrieved for each frequency when a cell flows in the channel. The shape of this pulse depends on the position and relative impedance of the microparticle and fluid in the channel, while its measured module and phase depends on the www.nature.com/scientificreports/ volume (or size) of the particle compared to the dimensions of the channel. When used as a spectroscopy with multiple excitation frequencies, it is possible to adequately characterize the dielectric properties of the cells flowing inside a microchannel, which can be used in association with their impedance models to retrieve important characteristics of microorganisms 28,29,45,49 . The impedance is generally converted to a complex permittivity since the geometrical parameters of the microchannel and electrodes are known, and used for the characterization. This is a simple and effective technique to count and characterize particles in a fluid, providing information in real-time for feedback control, or data for later analysis or post-processing 50 . The high volume factor between the cell and detection area when using coplanar electrodes in IFC creates a fringing effect between the electrodes, which is difficult to modelize. A simple empirical equation is given in 51 and 52 as an alternative to estimate the size of microparticles when using IFC. The particle diameter D can be estimated using a fit to the cubic root of the combined measured impedance magnitude difference | Z 1 | and | Z 2 | observed when a particle passes the first and second electrode pairs in the channel, respectively. This is shown in Eq. 2, where G is a constant that accounts for all the parameters linked to the electronics and fluidics, such as the electrode configurations, magnitude and frequency of the excitation signal, filter bandwidth, channel depth and width, electrodes width, EDL capacitance C edl , buffer conductivity, and electronics gains. The constant G can be determined empirically by testing the IFC system with circular beads of know diameters, then adjusting G until the estimated diameters match the effective beads diameters used during the test.
From the Randles 53 model and from the impedance model of single-shelled cells 42,43 , it is possible to deduce the optimal excitation frequency range for microorganism characterization, which is found to be between 100 kHz and 10 MHz. For frequencies lower than 100 kHz, the sensibility of the sensor to microparticles is reduced considering that the electrical double-layer (EDL) and ionic diffusion from the Warburg element dominate the measured impedance 38,44,46 . Above 10 MHz, the PCB dielectric begins to shunt the channel impedance, and the parasitics of the measurement electronics significantly reduces the precision of the results. Above 1 GHz, the dionic reorientation of water molecules also affect the measured impedance. Calibration algorithms can be used to compensate for the errors obtained when using excitation frequencies outside of this range.

Material and methods
The design and fabrication of the impedance-sensing system and microfluidic system will be described in this section based on the principles of IFC. Impedance-sensing system. The bloc-diagram of the impedance-sensing device is shown in Fig. 2. It is based on a LIA topology, and extracts the amplitude and phase of a high frequency input signal. It features two channels used to perform a differential analysis. A square waveform with frequency ranging from 200 kHz to 200 MHz is created by the clock-system (Si5351, Silicon Lab). This signal is sent to a quadrature generator to create two 90-degree phase-shifted square waveforms at half the clock signal frequency. The in-phase waveform is attenuated to 100 mVpp to keep a safe linear current response that will not harm the cells during the experiment, and then sent to the two differential electrode pairs of the microfluidics system. Two current responses are obtained, which are then amplified and converted to voltage signals by transimpedance amplifiers (TIAs). The www.nature.com/scientificreports/ outputs of the TIAs are mixed by two phase-sensitive detectors (PSDs). This yields four output signals representing the real and imaginary current responses of two electrode pairs, which can be sampled by the ADCs. The impedance magnitude and phase differences can be retrieved from the real and imaginary impedance components, and are processed by a square-to-sine spectroscopy algorithm to accurately retrieve the cells properties at a high throughput. The differential design is important to decrease the effects of the common noise and increase the sensitivity to the flowing particles 46 .
Square input signal. IFC and EIS systems are generally performed using sinewaves. The main advantage of square-waves over sinewaves is that they replace the complex hardware associated with conventional LIAs with much simpler clock-based circuitry 54 . For instance, the digital-to-analog converter (DAC), wave generator, and linear mixer can all be replaced by a simple clock-system and inexpensive controlled switches. This leads to a decrease of the system's power-consumption and cost, which are of a prime importance in portable applications. The square wave input voltage V in outputted by the clock-system and sent to the electrodes is defined in the time domain by Eq. (3). Square signals are multi-frequency signals. Looking at their Fourier transform, we see that a square wave is composed of a fundamental frequency followed by odd-frequency harmonics of decreasing amplitudes. An ideal symmetrical square wave of amplitude 2V o with peaks at −V o and +V o follows the geometrical sum of Eq. (3).
An important point to note from Eq. (3) is that the amplitude of the fundamental differs from the one of a sinewave of amplitude 2V o by a factor of 4 π . Such a difference biases the measurements and affects the precision but can be mostly corrected (up to a couple of percents) by following an algorithm proposed by Subhan 54 . Another consideration is the introduction of harmonics in the circuit, which raises the noise floor of the system.
Transimpedance amplifiers. Current responses are obtained from the two electrode pairs that follow Ohm's law for complex impedance. Current responses being difficult to interact with, two transimpedance amplifiers (TIAs) are used to convert them to voltage signals. TIAs are current to voltage converters generally implemented using one operational amplifier, as shown in Fig. 3. The practical implementation also uses a capacitor for stability in parallel with the resistor in the feedback loop. In the simplest case for a square input signal, Eq. (4) represents the voltage output of the TIAs in time V TIA (t) , where |Z| and θ are respectively the impedance magnitude and phase of the SUT for any harmonic angular frequency nω 0 .
To accomodate for a wide range of input impedances, a programmable gain array (PGA) with a feedback resistor R f and capacitor C f is added to the TIAs to control the gain at will. The PGA is achieved using a multiplexer toggled by a microcontroler that can switch between different gain resistors in the feedback loop of the TIAs 30 . The feedback capacitor is needed by the TIAs to prevent high frequency ringing. This can cause a limitation for high frequency measurements since an attenuation is expected at frequencies higher than a couple of megahertz because of the time-constant of the RC network formed by R f and C f . However, the prototype can still be used at higher frequencies with an adequate calibration, although with a reduced accuracy associated with the lessened measured signal amplitude. The trade-offs associated with TIAs are described in Orozco 55 www.nature.com/scientificreports/ Phase-sensitive detectors. The TIAs outputs cannot be sampled directly using an ADC since the frequency of the signals of interest is too high (the relevant harmonics can go as high as 110 MHz considering the five first harmonics of a 10 MHz square signal). A solution consists of using a mixing and filtering stage implemented from a phase-sensitive detector (PSDs). PSDs act as narrowband filter similarly to LIAs to precisely retrieve the amplitude and phase of a signal buried in noise 41 . PSDs use square signals and an inverter to switch between the original and inverted version of the signal of interest at the frequency of the square reference signal. This switching yields DC components proportional to the real and imaginary current of the SUT's impedance. The behavior and implementation of the PSD is shown in Fig. 3. The DC values of the real and imaginary components of the current responses at the output of the PSDs are described by Eqs. (5) and (6).
Quadrature generator. Operating the PSDs mixers requires two square signals in quadrature. Those signals can be precisely obtained from a quadrature generator circuit using a comparator and two D-Flip-Flops, as shown in Fig. 3. This technique is ultrawideband and relatively simple to implement but can be used only for low-power binary signals since the current is sunk directly from the low-power flip-flops. Programmable delays are added in the path of the reference signals to compensate for the delays of the TIAs circuits. This way, the measured phase response from the PSDs is only affected by the SUT.
Square to sine spectroscopy algorithm. Now, as can be seen from Eqs. (5) and (6), it is not trivial to recover the impedance magnitude and phase of the fundamental when using PSDs, as is the case with LIAs. Indeed, harmonics of the square excitation signal are present at every odd frequency of the fundamental which adds a systematic error to the impedance measurement 54,56 . The harmonics present in the square signal are multiplied together by the mixer and pushed to DC along with the desired fundamental frequency. This systematic error is non-trivial as it depends on the impedance response of the SUT. An algorithm inspired from Subhan 54 can be used to cancel the systematic error. The values of square impedance at the harmonics can be subtracted or added to the fundamental impedance following a certain set of rules www.nature.com/scientificreports/ described in Subhan 54 . It is thus necessary to measure the entire impedance spectrum before computing the corrected impedance at a given frequency. The real component of the impedance devoid of the systematic error V sine−ω 0 , follows Eq. (7) where E is the residual error after correction, which depends on the number of frequency points that were subtracted. There is, however, a practical limit to the number of points that can be subtracted considering that the impedance at that frequency must be measured (or extrapolated) beforehand, which might not be possible for high frequency samples. A similar process can be repeated for the quadrature component in Eq. (8). The corrected impedance magnitude and phase of Eqs. (9) and (10) can then be reconstructed.
Printed circuit board. Considering that the microparticles that pass in the channel are microscopic, the sensor has to possess a high sensitivity. For the electronics, a thorough understanding of noise and best PCB design practices is required. The impedance-sensing system is made from a four-layer PCB, and has a size of 50 × 50 × 15 mm including the components. The substrate is FR-4 TG150, with minimum spacing of 0.1524 mm and a thickness of 1.5 mm. Finally, the surface finish is HASL with 1 oz copper. The final PCB with all components is shown in Fig. 2.
Microcontroller unit. The IFC system uses a MSP430F5529 as microcontroller unit (MCU). The MSP430F5529 is a mixed signal MCU used in low-power applications. It dissipates about 6 mW when active and 24 µW when in low-power mode. 6 channels of the 12-bit ADCs are used by the impedance-sensing system to measure the real and imaginary values of the outputs of the phase sensitive detectors ℜ(V PSD−ω 0 ) and ℑ(V PSD−ω 0 ) of the two electrode responses, as well as the 5V power-supply voltage V DD and the battery voltage. The I2C module is used to program the clock-system and modify the excitation signal frequency. The UART module is used with the external integrated chip FT232RL to transfer the data to a nearby computer using Bluetooth or USB 2.0. A couple of I/O pins are used to enable the power-supplies and status LEDs, modify the gains of the PGA on the fly, reset the phase of the measurement by enabling or disabling the flip-flops of the quadrature generator, and reset the MCU.
Microfluidics system. The microfluidics system created for this study encompasses the micro-electrodes designed on PCB and a PDMS microchannel squeezed hermetically between 3d-printed components.
Microchannel. The microfluidic system manufactured in this study has an inlet and an outlet, where the liquid respectively enters and exits the device, fitted to soft Tygon thermoplastic tubing 57 . The inlet tube is linked to a glass syringe connected to a precise motorized syringe pump by Cole-Parmer model CP-120 that compresses the syringe at a constant programmable rate. The SUT flows from the syringe to the tubes before entering the inlet. It then reaches the PDMS microchannel where it is sensed by the PCB microelectrodes. The liquid finally exits through the outlet tube, which is connected to a waste container. The microfluidics system is shown in Fig. 4. The whole fabrication process is described in Bouzid 38 . A mold is initially drawn on a CAD software such as Solidworks. The drawn model is sent as a .STL file to the CADworks3D software to be meshed. This new meshed model is used by the stereolithography 3d-printer CADworks3d H50-405 to print a 3D-mold using Master Mold for PDMS Device Photopolymer Resin TM . The resin is then rinsed with IPA (90%) or methyl hydrate for 5 min and blow dried using an air gun. The mold is then cured using UV light around 400 nm in a LED light curing box for 50 min. PDMS (SYLGARD TM 184 Silicone Elastomer Base) and a curing agent (SYLGARD TM 184 Silicone Elastomer Curing Agent) are mixed at a ratio of 10:1 and degassed by letting the solution rest for 60 min. The mixture is put in the mold and cooked on a hot plate for 50 min at 70 • , then dried overnight at 40 • C 58 . A scalpel is used to gently prick-off the channel from the mold. The surfaces of the PDMS microchannel are then exposed to plasma at 600 KPa for one min. Polyethylene glycol is immediately applied on the surface to keep its hydrophilicity for longer periods of times. After a 10 min wait time, the PDMS microchannel is cooked for 10 min at 130 • C on a hot plate 58 . A dry gun is blown on the channel to take away any residues. The PCB electrodes are aligned on the microchannel and a pressure is applied to seal them. The PDMS channel and PCB electrodes are sealed tight by a system of 3d-printed squeezers that compress the channel and PCB electrodes. Those are tightened together by bolts, hermetically sealing the microchannel and PCB electrode due to the flexible nature of PDMS. Finally, tubes are inserted at the inlet and outlet of the PDMS microchannel, concluding on the whole process. The microfluidics system is thus hermetic and easy to handle. cos(θ(nω 0 )) n 2 |Z|(nω 0 ) = cos(θ(ω 0 )) www.nature.com/scientificreports/ The shape of the mold is shown in Fig. 4. The volume fraction (i.e. the ratio between the volume of the microparticle and volume of the fluid affected by the electrical field of the electrodes) should be maximized in IFC applications to obtain the highest sensitivity. Maximizing the volume fraction requires that the channel and electrodes size be identical to the microparticle of interest. In practice, the channel must be larger by a given safety factor to allow the fluid to circulate without clogging. Hemispherical bubble traps based on the work of Kang 59 were added to the channel to reduce the quantity of air bubbles reaching the electrodes which could falsify the measured impedance. Using the stereolithography 3d-printer H50-405, a theoretical resolution of 30 µm is possible, although in practice, the minimum size of a printable channel without major defects is 90 µm. For the case of this study, a channel size of 180 µm was chosen. The PDMS cured in this mold solidifies into a structure with four openings where are placed the bolts. The electrode and squeezers share the same shape and openings for the bolts as the PDMS.
Electrodes. Coplanar electrodes are chosen in this design because they are 2D structures that can easily be made using a lithography mask or directly on a PCB 28 . The non-homogeneous electrical field distribution of coplanar electrodes does introduce errors in the measurement since the cell's vertical position in the channel are subject to varying levels of electrical field. A higher particle in the channel typically experiences weaker electrical field than a low particle, which results in a lower perceived amplitude 51,52 , as was shown in Fig. 1. Since the amplitude of the spike is used to infer the cell properties, a significant error is thus observed. Three solutions can be used to counter this problem: (1) Using parallel facing electrodes placed diagonally opposed in the channel instead of coplanar ones 52 . (2) Using centering techniques such as dielectrophoresis, acoustophoresis, inertial focusing and sheath flows 51 . (3) Using coplanar electrodes with distinctive geometry to obtain additional information about the vertical position of the particle in the channel 41,51 . The 5-electrode configuration described in De Ninno 51 is one of those distinctive geometry and has been chosen for this study, and is shown in Fig. 2. The relative prominence of the signal obtained from such a configuration can be used to correct the measured particle size. The downside of such a technique is that it reduces the sensibility of the sensor since the intricate electrode geometry increases the sensing volume 60 .
The electrodes in this study are fabricated on a one-layer PCB. They have a size of 46 × 21 × 1.6 mm. The PCB uses the conventional substrate FR-4 TG130 in the exact shape of the microfluidics channel. The employed surface finish is immersion gold (ENIG) (1U") with 1 oz copper. Inert metals such as gold or platinum are used for the electrodes because of their convenience in casting for small dimensions, for their unlimited lifetime, and since other types of electrodes such as Ag/AgCl are unsuitable for high excitation frequency 28

Results and discussion
To test and calibrate the performances of the device, EIS analysis of discrete resistors and capacitors forming complex circuits were measured. EIS analysis on saline water, as well as the detection and characterization results of microbeads and buccal cells using IFC also follow. Finally, the performance of the system is discussed.
Complex impedance circuit. The performances of the impedance-sensing system are measured from the EIS analysis of a 10 k discrete resistor in series with a parallel combination of a 4.7 k resistor and a 100 pF capacitor. The impedance-sensing system has a lowest excitation frequency of 20 kHz and a highest of 12 MHz. The square excitation signal is initialized at the lowest frequency, samples 64 data points, then the frequency is incremented logarithmically until the end frequency is reached. When that is the case, the frequency is reinitialized to the lowest frequency, and the process begins anew. The impedance magnitude and phase are shown in Fig. 5, and were recorded for about 34 s at a sampling rate of 655 Sps, which amounts to about 320 data points per frequency. A lower accuracy and precision can be observed for low and high frequency. This is to be expected since the frequency range of the impedance-sensing system is defined only for frequencies above 100 kHz due to the LPF of the PSDs having a −3 dB low-pass cutoff frequency of 20 kHz, which does not totally attenuate the square excitation signal up until at least 70 kHz. The high frequency bias originates from multiple sources. Firstly, the op-amps used in the electronics of the numerous modules are limited in frequency by their slew-rates, input capacitance, and gain-bandwidth. Secondly, capacitive coupling is introduced because of the PCB dielectric, traces, and wires, which tend to attenuate the high-frequency components. Since square excitation signals are used, the harmonics at higher frequencies get attenuated until the signal behaves much more like a sinewave around 20 MHz. As a result, the perceived impedance from the impedance-sensing device is increased since the measured current response is attenuated at high frequency.

Figure 5.
Bode plot of the impedance magnitude and phase response of a 10-k discrete resistor in series with the parallel combination of a 4.47-k resistor and a 100-pF capacitor. 320 samples were taken from the SUT for each frequency, and the average and standard deviations are calculated and displayed on the error bar on the left. Four sets of data are displayed, the measured raw impedance, the raw impedance after calibration, the calibrated impedance after transformation using the square to sine spectroscopy algorithm, and the theoretical impedance of the SUT. www.nature.com/scientificreports/ To solve some of these issues and linearize the sensor, a calibration using a look-up table is realized using resistors of known values. Since the parasitics are singular to the electronics, the same nonlinearity will be found for different values of resistance, which can be used as a frequency-dependent factor to linearize the magnitude and phase curves. The square-to-sine conversion adapted from Subhan 54 is then performed on the calibrated dataset. The raw, calibrated and converted-to-sine results are shown in Fig. 5 for the median impedance at each frequency points. At high frequencies, a bias is observed both in the impedance and phase since higher frequency data points are not available to perform the square to sine conversion from Eqs. (7) and (8). A way to solve this issue would be to extrapolate the behavior of the system from the previous points and use that extrapolation in the square to sine conversion. For the sake of simplicity, no such correction was attempted in this study. Apart from that bias, errors of less than 3% are observed for the magnitude and phase, for the frequency range considered in the spectroscopy. This is comparable to the commercial devices presented in Table 1, at a fraction of the power consumption and cost.
Saline solution. Following the proof of functioning and the calibration, we measure the impedance spectrum of a complex system. A solution of saline water at 22 • C is passed in the 180 µm wide PDMS microchannel of the microfluidic system, and the EIS is measured. The calibrated and converted EIS curves for the two pairs of electrodes follow the behavior of a series capacitor and resistor, as expected from the Randles model. This dataset will be used to calculate the corrected impedance after using the square-to-sine transformation algorithm after the addition of microbeads and buccal cells.
The measured impedance varies slightly according to the pressure exerted by the liquid flow. This difference is caused by a slight contraction of the PDMS walls caused by the liquid pressure, which also increases the liquid volume measured by the electrodes. The same effect can be observed for variations in temperature of the liquid. Thus, the liquid pressure and temperature are controlled for the whole duration of the experimentation.

Microbeads.
In order to replicate more accurately the expected behavior from cells and microparticles, polyethylene microbeads are added to the previous saline solution. This SUT is kept at the same conditions as before, at 22 • C and passed in a 180 µm wide microchannel. Only one square excitation signal is used this time, at a frequency of 1 MHz. The corrected impedance can then be calculated with Eq. (7) using the measured impedance of the microbeads as fundamental and the EIS of saline water as the harmonics. Considering the small difference in impedance spectroscopy between both tests, errors less than 1% are expected.
To avoid overloading the MCU, the impedance is sampled at a fixed high-frequency rate of 5461 Sps, while data is saved to memory only when a significant difference is observed in the real parts of the measured impedance of either electrodes. When that condition is detected, a burst of 64 consecutive measurement points is saved. This method produces regularly fixed data point with dense bursts of data when an event is detected. This event detection can be caused by a microparticle passing in the channel, or by sudden changes in liquid property or microchannel geometry. Signal processing is performed offline to retrieve only the events associated with a particle detection. Firstly, the average of the magnitude difference is removed using wavelet decomposition, the signal is then de-noised, low-pass filtered, and smoothed so that the impedance spikes caused by the particles are easier to recognize. A particle detection algorithm is used on this dataset to recover the positions of the peaks. A peak detection algorithm is first used, followed by a weak supervision approach using Snorkel 61 to discriminate between the peaks obtained from microbeads, bubbles, or any other outliers. Most of the oddly shaped, or weird behaving particles are thus removed from the dataset automatically. With the particles peak locations, it is possible to recover the amplitude and width of the patterns, which are used to estimate the microbeads size using Eq. (2).
As an example, the pattern in (a) of Fig. 6 is studied. The first electrode has an impedance magnitude difference of 230 , while the second electrode has an impedance magnitude difference of 250 . Its G constant is estimated from the dataset to be around 10. This leads us to a diameter value of 78 µm. It is also possible to measure the particle velocity by dividing the distance between the electrode pairs L with the time it took for the particle to go from one electrode to the other t (which is the time difference between the two impedance maximums). The time it took for the particle to pass the electrodes is found to be 2.9 ms, while the distance between the electrode pairs is of 406 µm. This leads to a velocity around 14 cm/s. This flow rate has been found empirically to provide good and reliable measurements, since lower flow-rate can cause particles to stick to the walls of the channels or electrodes, and higher flow-rate are associated with a decreased time resolution. The fastest useful flow-rate for this sampling rate is when a minimum of 5 points are detected for a full particle. Any less than that is considered an outlier by the classifier. This leads to a maximum theoretical flow-rate of 89 cm/s. Buccal cells. The proposed system works to detect cells of a maximum size fixed by the width of the microchannel, and of a minimum size fixed by the sensibility and inherent noise of the sensor. For the case of this study, this leads to a minimum and maximum cell size of 45 µm and 180 µm respectively. An easy to test cell that fits those size requirements are those found in the mouths, the so-called buccal cells, with sizes typically ranging between 50 and 60 µm 62 . Those cells were scraped from the tongue and cheek of the corresponding author, and mixed with the same saline solution used in the previous tests. An example of a cell detection is shown in Fig. 6, where the 1st and 2nd pair of electrodes each detected successive impedance events, which were used to characterize the cell size around 51 µm.
Collected datasets. The true potential of IFC sensors lies in how automatizable the sampling and testing process can be for biological studies. It could be imagined that a team of biologists collects the impedance data of millions of cells and particles using the portable device described in this study, and efficiently extract the www.nature.com/scientificreports/ important information of these cells using machine learning and high-end signal processing. As a small-scale example, two datasets were collected from the proposed system and passed into the peak detection and classification algorithm. 447 beads were detected from 617 detection events, and 360 buccal cells were successfully detected from 2823 detection events. The debris in the saline solutions from the cheek and tongue scraping, and the low impedance difference obtained from the buccal cells explain why so much detection events were detected by the algorithm compared to the number of actual buccal cells. The diameters of the microbeads and buccal cells were estimated as described previously and compiled in the histograms of Fig. 6. A minimum impedance difference threshold is used to classify what counts as a detection event from the measurement noise. This lower threshold means that the microparticles of sizes below 45 µm that are present in the solution are not registered by the algorithms. It could be said that the effective sensitivity of the impedance-sensing system to detect small particles is of 45 µm when used with a 180 µm wide microchannel.

Measured performance.
The performance and characteristics of the presented impedance-sensing device and microfluidic system are summarized in Table 2. The impedance-sensing system created for this study is the first found in the scientific literature to achieve great sensitivity level over wide frequency and impedance range while boasting a small size, low-cost, and low power-consumption. The impedance-sensing device coupled with the microfluidics systems are effectively capable of measuring and estimating the properties of the microparticles of sizes going as low as 45 µm when used within a 180 µm wide microchannel. The dimensions of the microchannel are fixed by the limitations of the 3D-printer, which could be improved for this study by using a 3D-printer such as the one designed by Gong 63 . This homemade 3D-printer is specifically made for microfluidics and can attain truly microscopic scales of 18 × 20 µm by modifying the type of resin used and optimizing the stereolithographic process. This higher resolution would help increase the sensitivity of our device for smaller particle detection. The impedancesensing device takes 50 mm × 50 mm × 15 mm of space, while the microfluidics system is 46 mm × 25 mm × 50 mm, with a combined weight of 300 g, making them portable enough to be put in a backpack for applications in the field. The electrode pairs in the microchannel are separated by 424 µm and each have a width of 106 µm compared to the microchannel size of 180 µm × 180 µm. The impedance-sensing system only needs 1.2 W to  www.nature.com/scientificreports/ function adequately, and is powered by a low-voltage battery of 2.5-3 V. The power consumption of the system is sufficient for portable applications and could be powered for a couple of hours at a time. The power consumption could, however, be greatly reduced by creating a custom ASIC instead of using discrete components. Indeed, the vast majority of the power (about 80%) is dissipated in the op-amps, while they serve only to do basic functions such as inverting and amplifying signals that could be replaced by optimized high frequency transistors. The impedance-sensing system costs around $300, while the microfluidics system costs only $10 per microchannel excluding the initial cost of the 3D printer. The impedance measurement range between 200 and 120 k is similar to the portable impedance analyzer described in the literature, such as the ones from Al-Ali 37 , and Radil 34 . The frequency range is adequate for IFC applications, with the important frequency range between 100 kHz and 10 MHz covered by the device. The upper frequency limit of 12 MHz observed in this work is fixed by the limitations of the op-amps used in the TIAs. The limited bandwidth of the op-amps attenuates the harmonics of the square signal, which progressively modifies the square excitation signal into a sinusoidal shape. This introduces significant disparity for frequencies higher than 12 MHz which goes above the 3% precision reported for the device. The device can theoretically be used with excitation frequencies as high as 100 MHz, but the reported error would increase significantly. Finally, an excitation voltage of only 100 mVpp is used, which is low enough to not affect most microorganisms in that size range.

Conclusion
This study succeeded in creating an autonomous device for the characterization of microorganisms in the fields.
Using an inexpensive 3D printing manufacturing technique and standard printed circuit board technology, the presented device can detect and characterize microorganisms larger than 45 µm. The device succeeded in characterizing and differentiating between buccal cells and polyethylene microbeads. Future work will focus on improving the sensibility of the sensor to characterize microparticles of smaller sizes, as well as increasing the number of parameters that can be monitored to achieve a better characterization. Following the recent advances in micro-optical systems, adding a low-power 3D-imaging system to the device will be investigated.

Data availability
The Python, MATLAB, and C source code for this project, and the dataset obtained from the device are available upon reasonable request to the corresponding author. Python was used to sample the dataset from the device. C code was used to interact with the sensors and MCU. Post processing was done in Matlab and Python.